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We study the magnetic excitation spectrum of the spin-1 chain with Hamiltonian Tt = . cos 8Si ■ 
S i+ i +sin6>(Si ■ S 4+ i) 2 . We focus on the range < 9 < +7r/4 where the spin chain is in the gapped 
Haldane phase. The excitation spectrum and static structure factor is studied using direct Lanczos 
diagonalization of small systems and density-matrix renormalization group techniques combined 
with the single-mode approximation. The magnon dispersion has a minimum at k = tv until a 
critical value 9 C = 0.38 is reached at which the curvature (velocity) vanishes. Beyond this point, 
which is distinct from the VBS point and the Lifshitz point, the minimum lies at an incommensurate 
value that goes smoothly to k = 27r/3 when 9 approaches tt/4, the Lai-Sutherland point. The mode 
remains isolated from the other states : there is no evidence of spinon deconfinement before the 
point 9 = +7r/4. These findings explain recent observation of the behavior of the magnetization 
curve M « (# - H c ) 1/A for 9 = 9 C . 

PACS: 75.10.Jm, 75.20.Hr 



I. INTRODUCTION 



It is now well known that one-dimensional spin-5 Heisenberg antiferromagnets (AF) have qualitatively different 
properties according to whether the spin value 5 is integer or half-integeiu. The existence of a singlet-triplet gap just 
above the ground state is clearly establishedmJ by numerical techniques for the S=l and S=2 cases. Ot*r physical 
understanding of these phenomena is based on the construction due to Affleck, Kennedy, Lieb and TasakiQ (AKLT). 
These authors were able to obtain explicitly the ground state of the following bilinear-biquadratic Hamiltonian : 

"Haklt = ^ Si ■ Sj+l + -(Si • Si + i) 2 , (1.1) 

i 

where Sj are quantum 5=1 spin operators living on a chain whose sites are indexed by i. To construct the ground 
state wavefunction of the Hamiltonian (1.1), each original 5 = 1 spin is written as two 5=1/2 spins in a triplet 
state. Then the ground state is obtained by coupling into a singlet state all nearest- neighbor spins- 1/2, thus forming 
a crystalline pattern of valence bonds. This state is called the valence-bond-solid state (VBS). It is in fact the unique 
ground-state and excitations have a gap which is known rigorously to be nonzero. Many results followed from the 
VBS picture. Eos example, it implies that there are free spins 1 /2 at the end of an open chainQ. This has been verified 
experimentally&u by electron spin resonance on copper spins randomly introduced in the S=l AF chain compound 
NENP. The bulk excitations are easily pictured : by breaking a singlet bond into a triploUcoie creates a local objet 
that move along the chain. This is only an approximate eigenstate but it has good overlaalH~EJ with the first excited 
state which is a triplet, the "magnon" , as predicted in the original derivation of Haldanetl. 



This appealing picture is certainly correct for the Hamiltonian (1.1) but it remains to be understood how closely 
it applies to the standard Heisenberg exchange Hamiltonian for which the biquadratic coupling in Eq.([0]) is zero. 
Initial studieso proposed the generalized family of models : 

Hg = J2 cosdS * ■ S *+i + sintf (S, • S 4+1 ) 2 . (1.2) 

i 

This family includes the familiar Heisenberg model for 9 = and the AKLT Hamiltonian for tan6Vss = 1/3. Since 
the gap does not go to zero in the interval [0, Ovbs] it is likely that there is no phase transition and that the two 
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limiting Hamiltonians 8 = and 8 = 8vbs — 0.3218 share the same physics. It is exactly the same line of arguments 
that justify— the use of Laughlin wavefunctions to describe correlated states of electrons in the Fractional Quantum 
Hall EffectH 

However, the AKLT point has special properties : for example the spin correlation functions have a spatial decay 
which is purely exponential, co ntra ry to the Ornstein-Zernicke decay which is observedJbr the Heisenberg model. A 
study of the family of models (1.2) revealed that the AKLT point is a disorder pohrO beyond which short-range 
incommensurability appears in real space spin correlationsEZl Hence, it is at a boundary in some sense. The evolution 
of the fPffn excitations as a function of 8 is the subject of the present paper. For 9 = ir/A it is known from the Bethe 
AnsatzcSEj that the spectrum is gapless and that excitations form a continuum as in the case of the spin- 1/2 AF 
chain. Near the Heisenberg point, it is known that there is an isolated branch of triplet excitationstj that enter in 
a two-particle continuum for a wavevector q sa 0.37T. In the neighborhood of q = 7r the continuum is well above the 
isolated mode. This holds also for the AKLT Hamiltonian, as shown by Lanczos diagonalizatiorJ 2 -! Since there is no 
phase transition with zero gap between 8 = and the Lai-Sutherland point 8 = +7r/4, the common belief is that 
things evolve smoothly. This rather vague statement deserves however more scrutiny in view of the change of the spin 
correlations at the AKLT pointO. An interesting issue is the fate of the Haldane triplet mode that should ultimately 
disappear in the spinon continuum for 8 = tt/4. It has been speculated recentlyea that beyond the AKLT point 
8 > 8vbs the Haldane mode disappears and is replaced instead by a gapped spinon continuum which becomes gapless 
only right at the Lai-Sutherland point. However this intriguing picture is baaed only on variational wavefunctions 
whose relevance to the problem is at least unclear. It has also been observed^ that the magnetization curve M(H) 
displays intriguing behavior for some value of 8. Some recent workEa has undertaken the study of dynamical properties 
between the Heisenberg and VBS points. 

In this paper, we investigate the excitation spectrum for a range of values of 8 in the Haldane phase. Our main 
finding is that the Haldane mode remains well-defined and isolated from other states : there is no deconfinement of 
spinons. The minimum wavevector of the dispersion relation E(k\ stays at the zone boundary k = n till a critical 
value 8 C w 0.38 which is between the 9vbs and Lifshitz pointsE3 8l and clearly different from these two values. 
Beyond 8 C the dispersion has a minimum at some incommensurate wavevector that evolves smoothly towards 2-7r/3 
when 8 — > 7r/4. At 9 C , the dispersion curve has a quartic minimum and we show that this implies a critical behavior 
w (H - H C W 4 for the magnetization M(H). 

In sect. (O), we describe our results for the dispersion of the Haldane triplet obtained from Lanczos exact diagonal- 
izations. In sect. (HI), the single-mode approximation is applied to DMRG calculations of the static structure factor in 
order to estimate the lower edge of the excitation spectrum. We discuss the physical properties of the Haldane phase 
and the relationship with the nonlinear sigma model in sect. (IV). We discuss also the shape of the magnetization 
curve in sect. (IV). Finally our conclusions are presented in sect. (M) 



II. LANCZOS RESULTS 



At the Heisenberg point 8 = 0, the excitation spectrum of the hamiltonian (1.2) has been studied in considerable 
detail. We first use the Lanczos diagonalization method to study the interval < 8 < +ir/4. We use conservation 
of the z-component of the total spin and compute the eigenenergies of the first 10 low-lying states. We use also 
translation symmetry to work in sectors with fixed lattice momentum k. To keep the iiamiltonian real even when 
the momentum k is not or ir, we use a non-trivial change of basis due to TakahashiES. This requires up to 350 
Lanczos iterations which is considerably more expensive than just getting the ground state of each S z sector. To 
prevent runaway of the Lanczos process we are forced to re-orthogonalize with respect to the previous basis vector at 
each step. This bottleneck leads to a limit on the length of the chain which is 16 sites. This is a difficulty which does 
not appear when computing only the ground state. Our results for 8 = are shown in Figure [l]. Above the ground 
state, there is a well-defined triplet branch which is the prominent Haldane mode, the "magnon" . The. minimum of 
the dispersion is at the zone boundary k = n. This is the Haldane gap. The mode enters a continuumEj at k 0.37T. 
At k = the continuum is the bottom of the spectrum and starts at twice the Haldane gap (approximately for N=16 
sites on our figure but checked with excellent precision). 

There is a nonzero curvature of the dispersion E(k) at k = it. We model the bottom of the dispersion relation by 
a rclativistic formula as suggested by the nonlinear sigma model : 

E{k) = v/A 2 +w 2 (fc-fc ) 2 . (2.1) 

Here ko is wavevector of the minimum, A is the gap and v is the velocity. Close to the minimum we have : 

£(fc)«A + ^(fc-fc ) 2 - (2.2) 
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The second derivative of the dispersion is thus simply related to the velocity that occurs in the nonlinear sigma model 
description. 

When we increase 9, there is a range of values for which nothing qualitatively new happens. This range includes 
the VBS point 9 = Qvbs — 0.3218. The spectrum at this point is displayed in Figure ^. It shows the same features as 
that of Figure [l| Notably there is still a nonzero curvature at k — tt of the magnon branch. If we consider the energy 
of the state at fco — t an d the two closest states on the magnon branch with k = tt — 2ir/N and k = tt — 4ir/N, we 
can fit by a fourth-order polynomial : 



c • \2 , d n. i. , i 



E(k) = E(k ) + -(fc - koY + — {k - fc ) 4 . (2.3) 

and hence obtain an estimate of the velocity. Close to the VBS point, the values of c as a function of the chain 
length are given in Figure |^. At this point it is important to note that the finite-size effects are in fact extraordinarily 
small. This is related to the fact that the spin-spin correlations in the VBS wavefunction are of very short range 
£ = 1/ log 3 » 0.9. So our estimate for c is very precise : from N=4 to N=16 sites the dependence is very weak. We can 
exclude a vanishing velocity at this point. Indeed we find c = 0.9778(1). Similarly the Haldane mode stays isolated 
from the continuum which lies much above : these states display extremely fast convergence to the thermodynamic 
limit. 

We observe in Figure || that the velocity vanishes at a point 9 C = 0.38 which is between the VBS point 9vbs — 0.32 
and the Lifshitz point 9l = 0.414. There the minimum of the dispersion relation becomes of fourth order but still 
there are no states coming from above to fill the gap between the Haldane mode and the continuum. This case is 
shown in Figuxe ^. Once again we are in a regime of very weak finite-size effects : measurements of the spin correlation 
in this regimeEH are in agreement with the present findings. 



We can also obtain an estimate of the fourth derivative d of the dispersion relation E(k) through Eq. (2.3). It is 
given in the neighborhood of the VBS point in Figure |^. The size dependence is minimal again at 9vbs- Note that 
right at 9 C the fourth derivative d clearly extrapolat e to a positive non-zero value : this fact will be important to 
explain the shape of the magnetization curve in sect.(|Tv|). 

Beyond the special point 9 C , the minimum of the dispersion no longer lies at fc = n but is at a wavevector which 
is in general incommensurate and dependent upon 9. The curvature of the dispersion is now negative at the zone 
boundary : see Figure |^. A typical case is shown in Figure ^ for 9 = 0.52. There is no evidence for the dcconfined 
spinons proposed by Yamamotocj. If this was the case, the corresponding energy states should be filling the empty 
interval above the Haldane triplet mode. We see no evidence for this. 

When 9 — > ir/4 the gap closes and the overall picture is consistent with the apparition of the tripled period 27r/3. 
However in this regime finite-size effects become important as the phase transition is approached. It is known that 
right at the Lai-Sutherland, point, the phase transition is described by a a SU(3) generalization of the Kosterlitz- 
Thouless phase transitionEZI For 9 > 7r/4, the system eaters a gapless phase. In the regime — 7r/4 < 9 < 7r/4, it is 
extremely likely that the gap vanishes only at 9 = ±7r/4lI3lI3. However, we cannot exclude from our numerical study 
that spinon deconfinement appears before that point but this would involve yet another phase transition for which 
there is at the present time no evidence. 



III. SINGLE-MODE APPROXIMATION 



To confirm the results obtained in the previous section, we now turn to the single-mode approximation of the 
dispersion relation. This approximation was introduced originally by Bijl and Feynman to determine the phonon- 
roton dispersion in superfluid 4 He. It has been also successfully applied to Haldane gap systemstLil. 

The dynamical structure factor of the spin system is defined by : 



S(q,w)= dte^^e*« n (SZ(t)SZ). (3.1) 

J-oo n 

In this equation, the a index need not be specified since we consider only isotropic systems. The equal-time correlation 
function is : 

S(q) = (S«S%) = A- f +0 ° duS(q,u>), (3.2) 

where = e lqn S*. Ima gine that we know the ground state |0). Then a guess for the first excited state 

may be simply |0). This is a triplet with lattice momentum q. This will be close to an exact excited state if the 
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dynamical structure factor S(q,uj) is strongly peaked at the energy wsma(?) of the state Sg\0). This energy is given 
by the formula : 

oJsuAiq) = 2^)<[ 5 -<> (3.3) 
where Tig is the Hamiltonian. The a dvan tage of the single-mode approximation is that it does not require a full 



dynamical study of the system : in Eq.(3.3), it is possible to evaluate all quantities once we have a good approximation 



for the ground state. Th is is feasible by the DMRG algorithm which is extremely efficient in quantum spin chains. 



The commutators in Eq^S^) can be evaluated straightforwardly with the result : 

u S MA{q) = (cosq - l)c(6)/S(q), (3.4) 

where c(9) is some constant given by : 
c*(6) = cose [< SfS* +1 > + < S?S? +1 >] + sin0[2 < (S. t ■ S l+1 ) 2 >-<($■ S i+1 )Sf + SfSf +1 (Si ■ S t+1 ) > 

+ 2 < s?s*s* +1 s? +1 + s*s?s? +1 s* +1 > 2 < (sff(s? +1 f + (svf(s* +1 f >]■ (3.5) 

We find that c(9) only varies slowly with 9 and it is now obvious that at least within the single-mode approximation 
the minimum in u>(q) will start to move away from tt before the Lifshitz point is reached (at the Lifshitz point the 
maximum of S(q) moves away from n). At the VBS point the ground state \VBS) is simple enough that ojsmaQ^) 
can be calculated exactlyS. One has uJsma 

(k) = J^L(5 + 3cosfc). Hence, we see explicitly that at the VBS point 
the second derivative and hence the velocity is non-zero within the SMA approximation. The value of the second 
derivative csma — 1-054 is very close to our Lanczos estimate of sect. II c = 0.9778. 

Using the above expressions, c{9) as well as S(q) have been calculated with the DMRG. We present results in the 
interval between 6 C as defined in sect. (jl^) and 9l- In Figure]]], the structure factor is plotted in the neighborhood 
of the zone boundary q = n. It becomes very flat, i.e. of fourth order at the Lifshitz point 9l — 0.414. However, 
the minimum in u>sMA{q) starts to move away from it for 9 S [0.38, 0.39] : this can be seen on Figure ^ where this 
quantity is plotted close to q — ir. This is consistent with our findings from the exact diagonalization results. Clearly, 
9 C and 9l are distinct, and in fact a considerable interval occurs between these two points. 

Inherently the SMA approximation assumes that 

S(k,w) = S{k)5(uj-E(k)) 

and hence neglects any incoherent background. Thus, LosMA(k) always over estimates the true excitations. This is 
important to remember when one compares the SMA results for the gap with other numerical estimates. Roughly 
the distance that the SMA value is above the true gap measures how big the incoherent contribution to S(k, ui) is. It 
works quite well around the Lifshitz point. 



IV. STRUCTURES WITHIN THE HALDANE PHASE AND MAGNETIZATION PROCESS 

From the numerical studies, both Lanczos and DMRG for single-mode approximation, we obtain a rather intricate 
structure within the Haldane phase — 7r/4 < 9 < +w/4. 

• In the interval — 7r/4 < 9 < 9vbSi the dispersion relation has its minimum at k = it, the spin correlations are 
commensurate : (So • S x ) ~ (— ) x exp(—x/^)/^/x. In this regime the 0(3) non-linear sigma model is a perfectly correct 
description of the low-energy long-wavelength behavior of the system in agreement with the original derivation by 
Haldanetl. 

• In the interval 9vbs < 9 < 9 c , the dispersion-Delation still has its minimum at k — ir but now the spin correlations 
oscillate with a period which is incommensurateEj, This is not seen in the static structure factor S(q) which remains 
peaked at q — it, a feature due to the short-range nature of the spin order. While E(k) is still sigma model like, this 
regime is not described by the sigma model which has only commensurate spin correlations. Indeed the nonlinear 
sigma model uses the Ansatz : 

S; = (-)*& + L i5 (4.1) 

where the fields <fr and L are smooth i.e. have Fourier modes only near k = 0. This is no longer valid even without a 
phase transition as a function of 9. 
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• In the interval 9 C < 9 < 9l, the dispersion has now a minimum away from the zone boundary. In fact there are 
two nonequivalent wavevectors 7r±Q at the minimum. The quantity S jq jl i s still peaked at q = tt. A possible effective 
theory would be now a non-linear sigma model describing helical orderc3Ea. It would involve a 3 x 3 rotation matrix 
as an order parameter to describe the short-range spin order. It is a phase with incommensurate short-range spin 
correlations and the spinon are confined. 

• In the interval 9l < 9 < +7r/4, the only difference with the previous case is that there is a double-peak structure 
in S(q). 

The above observations have interesting consequences fot,±he magnetization curve M(H). For the Heisenberg model, 
its shape and essential characteristics are well understoodEM When the applied uniform magnetic field H is less than 
the Haldane gap, there is no net magnetization : since H is coupled to the z-component of the total spin, the Zeeman 
Hamiltonian commutes with the Heisenberg Hamiltonian and hence does not change wavef unctions. So the ground 
state is unaffected by H. However at a critical value H c i = A, the Haldane gap, there is a crossing of levels and a 
magnetized state becomes the ground state. Finally there is full saturation to a ferromagnetic state beyond some 
field H C 2. Near H c i, the critical behavior is given by M{H) ~ [H — H c \) x l 2 . The exponent 1/2 is directly related 
to the dispersion relation of the magnons near k = n. Indeed the magnetization is given by inverting the relation 
H — H c \ = dE/dM and the function E(M) is the same as for a system of noninteracting fermions : 

£ = (A- ff )M + ^ (42, 

In this equation, N is the number of sites and v is the velocity which is defined through the dispersion relation : 

v 2 

E(k) = A + — (k - it) 2 + c 4 (fc - tt) 4 + 0(k - tt) 6 . (4.3) 

Since M = Nkp/n, one has : 

E(M) = M(A - H) + ^--^ + d 4 M 5 + 0(M 7 ). (4.4) 
6A L z 

Near the critical field this leads to M ~ (H — A) 1 / 2 . This depends crucially on the quadratic dispersion relation. If 
now we consider the generalized Ha milt onian at 9 = 9 C , one has v — but the fourth-order derivative is nonzero as 
can be seen in Figure ||. Hence Eq.(4.3) should be replaced by E(k) = A + c±(k — 7r) 4 . T here is no cubic term by 
imposing invariance under k^2n—k which is parity and lattice periodicity. Now in Eq.(4.4), the M 4 term dominates 
and this leads immediately to M ~ (H — A) 1 / 4 . Our findings explain the observation of such a behavior by Okunishi 
et alB3. When 9 > 9 C the dispersion has again a quadratic minimum and the line of arguments with v ^ is agaia 
valid : M ~ (H — A) 1 / 2 is the correct behavior in this range. The massless phase above H c \ however is likely to beU 
a two-component Luttinger liquid, contrary to the one-component Luttinger liquid that arises for 9 < 9 C . 



V. CONCLUSION 



We have studied the magnetic excitation spectrum of the bilinear-biquadratic S=l Hamiltonian that includes the 
Heisenberg point as well as the VBS point. We have shown by Lanczos and DMRG techniques combined with the 
single-mode approximation that the magnon dispersion becomes incommensurate at a critical value 9 C — 0.38 which is 
different from all previously known points in the phase diagram. Right at the special value 9 C the magnon dispersion 
has a fourth-order minimum at the zone boundary k = ~k. As a consequence the magnetization curve M(H) near the 
lower critical field is drastically modified : instead of a square-root behavior, the magnetization rises with a power-law 
1/4. This explains recent observations. 

From our study it is now clear that there are hidden structures within the bulk of the Haldane phase. This phase 
extends in the interval — 7r/4 < 9 < +7r/4 and the gap vanishes at the two boundary points. These critical points 
correspond, as is well known, to transitions towards completely different phases. It is widely believed that absence of 
phase transition means that the physics evolves smoothly, hence, that the Heisenberg point 9 = and the AKLT point 
9 = 0.32 share the same physical properties. It was known that spin correlations do change qualitatively even within 
the Haldane range — 7r/4 < 9 < +7r/4. We have shown that also the dynamical properties do change qualitatively 
and not only quantitatively. Curiously enough the special point characterizing the change 9 C is distinct from the 
VBS point. The bilinear-biquadratic S=l chain is thus an interesting counterexample showing the variety of physical 
phenomena that occur within a massive phase. 

Finally we have shown evidence that there is no spinon deconfinement when we approach the SU(3) phase transition 
at 9 = +7r/4. Thus, the Haldane phase is an example of a gapped spin liquid with incommensurate spin excitations, 
a phenomenon that may also appear in the normal phase of the underdoped cuprates. 
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FIG. 1. Excitation spectrum of the N=16 chain for the Heisenberg S=l AF chain, 9 — 0. Energies are relative to the ground 
state. The number denote the total spin value of each state. We have introduced small horizontal shifts to avoid superposition 
of symbols. The energies are given as a function of the chain momentum k (only the first 10 levels are calculated for each 
momentum and S z ). Note t he p rominent triplet branch that starts from the Haldane gap at the zone boundary at k — tt. The 
solid curve is the fit by Eq. (2.3). 
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FIG. 2. Excitation spectrum of the N=16 chain for 9 = Ovbs- The physics is qualitatively similar to 9 — 0, there is evidence 
for a nonzero curvature of the dispersion at k = it of the magnon and there is no spinon continuum. 
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2nd Derivative 
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FIG. 3. The second derivative c of the dispersion relation computed with the triplet state at k = n and its two closest 
neighbors along the magnon branch. Data are for several sizes N=8,10,12,14,16. It is nonzero at the VBS point but vanishes 
for 6 C = 0.38. Note that finite-size effects are extremely small in this range. We find cvbs = 0.9778(1). 
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FIG. 4. Same spectrum as figure 1 but for 9 
minimum at the zone boundary. 



6 C . At this point the dispersion of the triplet mode has a fourth-order 
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FIG. 6. Same spectrum as figure 1 for 6 = 0.52. This point lies deep in the incommensurate regime. The minimum wavevector 
is now inside the Brillouin zone and the magnon mode remains isolated : there is no evidence for spinon deconfinement. 
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FIG. 7. The static structure factor from a DMRG calculation. It is drawn only near the zone boundary for several 9 values 
in the interval 6 c ,6l- 
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FIG. 8. The dispersion LusMA(k) obtained from the single-mode approximation. The minimum moves away from k = 7r at 
the value 9 C which is the same as the Lanczos value 0.38(1). Energies have been divided by cos#. 
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